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We discuss how various models of scale-free complex networks approach their limiting properties 
when the size A'' of the network grows. We focus mainly on equilibrated networks and their finite- 
size degree distributions. Our results show that the position of the cutoff in the degree distribution, 
fccutoff , scales with A'' in a different way than predicted for N ^ oo, that is subleading corrections 
to the scaling fccutoff ~ are strong even for networks of order A'' ~ 10^ nodes. We observe 
also a logarithmic correction to the scaling for degenerated graphs with the degree distribution 
7r(fc) ^ k~^. On the other hand, the distribution of the maximal deg may have a different 

scaling than the cutoff and, moreover, it approaches the thermodynamic limit much faster. We 
argue that fcmax ~ N°' with the exponent a' = min(a, 1/(7 — 1)), where 7 is the exponent in the 
power law 7r(fc) ~ k~^ . We present also some results on the cutoff function and the distribution of 
the maximal degree in equilibrated networks. 



I. INTRODUCTION 

Recent progress in understanding the structure and 
function of complex networks [l| has been largely in- 
fluenced by the application of statistical methods of 
modern physics. The statistical mechanics of networks 
0, H, 0, H, B 0, Hi, if restricted to structural proper- 
ties, deals with two classes of problems. In the first 
one, one considers networks being in a sort of equilib- 
rium 0, 0, S, Si where the concept of statistical ensem- 
bles applies in a natural way. On the opposite, growing 
networks are usually treated by rate-equation formalism 
0, @|- Both approaches, however, share one common 
thing: the majority of considered models is analytically 
solvable only in the thermodynamic limit. This causes 
several difficulties. First of all, a question arises about the 
statistical equivalence of different ensembles of networks 
like canonical or grand-canonical ensemble Second, 
considering dynamical processes on networks, one can ask 
how they are infiuenced by approaching the thermody- 
namic limit in different ways [3|. 

The third question, which we want to address in this 
paper, is how fast is the convergence towards the limit- 
ing values of some network features, for different models. 
This question is much more important than the analo- 
gous question in more traditional physics dealing with 
atoms or molecules, where the number of degrees of free- 
dom is typically of order 10^^. Here, on the contrary, 
the largest observed system is the World-Wide Web with 
jV ~ 10^ — 10^° nodes. Many real networks are much 
smaller than this, typically having only lO'^ — 10* nodes. 
Therefore, finite-size corrections to solutions obtained in 
the thermodynamic limit are much stronger and in many 
situations cannot be neglected. 

As a most prominent example we will consider the de- 
gree distribution 7r(fc) which is of most importance for 
characterizing the network as well as many dynamical 
processes taking place on it. Finite-size effects are es- 
pecially strong for scale-free networks, for which 7r(fc) ~ 



k~'^ exhibits a power law. Examples are the Internet, the 
WWW and many social and biological networks. Since 
the power law cannot extend to infinity for finite size 
N, the degree distribution Tr{k) must decrease rapidly 
above some characteristic value k = fccutoff- The cutoff 
scales typically as fccutoff ~ N"', where the exponent a 
depends on the network type. For 2 < 7 < 3 when the 
second moment is divergent in the limit of large N, the 
cutoff appears either directly, or via moments of 7r(fc) in 
many applications, for instance in percolation or infec- 
tion spreading [l3|, transport models synchroniza- 
tion processes and others, see Ref. [l3|. Because qualita- 
tive differences may show up for various values of a, it is 
important to know the proper scaling of fccutoff with N. 

For a growing network like the Barabasi- Albert (BA) 
model [l3| or its generalization, the growing network with 
re-directio n ( GNR) model Q, the scaling is well estab- 
lished [l5| : for any 7 > 2 , the exponent a = 1/(7 — 1). 
Moreover, it can be shown that for TV large enough n{k) 
can be approximated as noc{k)w{k/N°'), where 7roo(fc) is 
the degree distribution for iV 00 and w{x) is some 
function independent of N. This approximation holds 
even for quite small networks of order lO'^ nodes because 
the convergence towards the limiting distribution is fast. 

The situation is not so clear for equilibrated net- 
works, that is networks in which evolution is governed 
by rewiring of existing connections rather than by adding 
new nodes, and thus can be regarded as being in a sort of 
equilibrium. Although the scaling exponent a has been 
estimated for some models [H, [3, [l^, the cutoff 
function w{x) has not been determined so rigorously as 
for growing networks. 

In this paper we will show that the predicted scaling 
is far from being true even for quite large networks. The 
paper is organized as follows. In Sec. II we discuss three 
models of equilibrated networks whose large- iV behavior 
we want to study. Section III starts with basic concepts of 
how to calculate finite-size degree distributions and how 
to extract the behavior of the cutoff function. Then we 
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consider two models, for which precise, semi-analytical 
results are available for sizes up to ~ 10^. In Sec. 
IV we discuss a more compHcated model and show how 
to simulate it on a computer. Section V is devoted to a 
relation between the cutoff and the maximal degree. The 
paper is closed with a short summary in Sec. VI. 



II. MODELS AND THEIR PROPERTIES IN 
THE THERMODYNAMIC LIMIT 

We shall start from defining three different models of 
equilibrated networks, whose finite-size properties will 
be further examined. The word "equihbrated" means 
that networks (graphs in mathematical language) are 
maximally random under given constraints. These con- 
straints are what defines the statistical ensemble of net- 
works. The statistical ensemble consists of a set of 
states - graphs {gi}, and a set of their statistical weights 
{Vr(ffi)}- This means that every graph g has a probabil- 
ity of occurrence proportional to W{g). Every physical 
quantity X is then defined to be the average over the set 
of graphs: (X) = E,, W^(3)^(.9)/ Eg ^^(.9)- Changing 
the set of graphs and/or the weights, one can obtain dif- 
ferent models of random graphs. Three such ensembles 
will be considered in this paper and are defined below. 
In all cases, we assume that the graphs are undirected, 
have labeled nodes, fixed number of nodes TV and links 
L, and that the weight W{g) of every labeled graph has 
a product form: 



N 



Wig) = n^'^^')' 



(1) 



where ki is the degree of the ith node of graph g and 
p{k) is some arbitrary node- weight function. The product 
measure H]) turns out to be very convenient to obtain the 
desired degree distributions in the thermodynamic limit 
by tuning the function p{k). In what follows, we will 
assume that all models have the same distribution 7r(fc) 
in the thermodynamic limit. Different models will be 
specified by restricting the set of possible graph shapes 
and/or p{k). 

Equilibrated simple graphs. We consider graphs 
without self- and multiple connections. Denoting by 
■Koo{k) the degree distribution for iV ^ cx3, we have: 



(2) 



as follows from Refs. [3,[23|- The constants A, B are cho- 
sen to have E/t T^oc{k) = 1 and Efe fci'oo(fc) = '2.L/N = k. 

ways. First, 
pik) and the 
calculate the 
degree distri- 
In the latter 
ye degree k is 



Equation ^ can be used in two different 
it tells us what 7roo(fc) would be for given 
average degree k. Second, it allows us to 
weight p{k) which will produce the desired 
bution TTdesik) in the thermodynamic limit, 
case, if N, L are chosen so that the averae 



equal to (fc) = Efe ^''^dcsik), then B = and the choice 
p(k) = 7rdos(fc)fc! implies that iToo(k) = TTdes(k). 

Equilibrated multigraphs. We again assume the 
product weight H]), but now we accept also degener- 
ated graphs, that is graphs with multiple- and self- 
connections. It can be shown 0,[l3;[2l| that the partition 
function of the system, being the sum over all configura- 
tions, assumes the form: 



z{n,l)=j2---y: 



p{ki) p{kN) 



fci=0 



fcjv=0 



Cat 



— 02L,fei + ---+fcN, 



(3) 

so it is equivalent to that of the balls-in-boxes model 
[131 or the zero-range-process model [l^l with weights 
p(k)/k\. The same formula ([2]) as for simple graphs holds 
for the degree distribution in the limit N ^ oo. 

Equilibrated trees. This ensemble consists of all la- 
beled, connected tree graphs with N nodes. In addition, 
we assume that one node is distinguished by a "stem" 
attached to it, which is convenient from a mathematical 
point of view, but it does not change the large-iV be- 
havior. Such "planted" trees can be treated in a special 
way and a number of quantities can be calculated analyt- 
ically [iqI. [23|. For instance, in the thermodynamic limit 
the degree distribution is given by 



pjk) 

{k-iy: 



(4) 



Note that we have (fc — 1)! in the denominator, in con- 
trast to Eq. (121). Now, the desired degree distribu- 
tion 7rdcs(fc) can be obtained only if Efe ^''''dcs(^) = 2, 
because the average degree for trees approaches 2 for 
N ^ oo. When this criterion is fulfilled, then assum- 
ing p{k) = 7rdcs(fc)(fc - 1)! one obtains 7roo(fc) = 7rdcs(fc)- 
All these models share one common property: the de- 
gree distribution 7roo(fc) in the thermodynamic limit is 
proportional to p{k). Therefore, by choosing p{k) ~ 
k\k~'' (for graphs) or ~ (fc — l)!fc~''' (for trees), one can 
make these networks scale-free. In this paper we shall 
stick to the following choice for the degree distribution 
in all models: 

(7 ~ l)r(27 - 3) r(fc + 7 - 3) 
''^^'')- r(7-2) r(A: + 27-3) ' 

for A; > 0, and 7roo(0) — 0. This is precisely the degree dis- 
tribution TToo (fc) in the GNR model of a growing tree men- 
tioned above. The average degree for this distribution is 
(fc) = 2, therefore the equilibrated trees will approach 
7roo(fc) in a natural way if p{k) = TToa{k){k — 1)!. As to 
equilibrated simple graphs and multigraphs, one then has 
to ensure that p{k) = 7roo(fc)fc! and k = 2L/N — > 2 for 
N oo, which can be simply done by assuming L — N. 
The purpose of choosing this particular distribution is 
that since finite-size effects in the GNR model are known 
we can compare what happens if networks are equi- 
librated but have the same 7roo(fc) as growing ones. 
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III. DEGREE DISTRIBUTION FOR A FINITE 
NETWORK 

In the previous section we discussed the behavior of the 
three models in the thermodynamic limit. Now we shall 
ask, how the degree distribution looks like for N < co. 
Assume that TrN{k) and 7roo(fc) are degree distributions 
for finite N and iV ^ oo, respectively. It is convenient 
to write the finite-size distribution TTjs[{k) as a product 
of 7roo(fc) being A^-independent, and some cutoff function 
w{N, k) depending explicitly on the size N: 



TTN{k) = nooik)w{N,k). 



(6) 



The function wlN, k) is model dependent. In Ref. jl^ 
it has been found that for the GNR model, moments of 
w{N, k) scale as follows: 



^u.(^,fc)fc" oc 7V"("+i)(l + 0(iV-")) (7) 



with the cutoff exponent a defined in the introduction 
and equal to a = 1/(7 — 1). This means that for suffi- 
ciently large TV the cutoff function depends effectively on 
a single rescaled variable x — k/N"": 



w{N,k) ^ w{k/N"). 



(8) 



For the GNR model, it is possible to find an expHcit form 
of w{x) for some values of 7. It is usually complicated, 
but it always has the following large- a: behavior: 



hiw{x) 



-X 



1/(1-") 



(9) 



where we defined the exponent 77 = 1/(1 — a). Formu- 
las lO, (HI) and ^ were verified numerically for N 
of order 10^ — 10"*, and an excellent agreement has been 
found. In this paper we ask, to what extent these rela- 
tions are valid for equilibrated networks. In particular, 
if the scaling ^ still holds, the moments will behave 
as /im ~ 7V"(™+i) for very large networks. The conver- 
gence towards this asymptotic behavior can be, however, 
different from that of growing networks. Therefore, we 
can relax the constraint that the sub-leading term de- 
cays with the power a and generally expect the following 
large- iV behavior: 



(10) 



with some constants a > 0, 5 7^ and /3 > 0, so that 
the subleading correction decays like with (3 not 

necessarily equal to a, as it was in Eq. 



A. Multigraphs 

We shall start from multigraphs. As we said, we 
assume that the average degree k is chosen to ensure 
J2k kiTooik) = k. This means that for a given number of 
nodes TV, the number of links L = L{N) is fixed. In case 



of the distribution L = N . The partition function 
^ becomes a function of TV only and can be rewritten 
as: 



^F-(z), 



(11) 



where the contour of integration encircles zero and F{z) 
denotes the generating function for 7roo(fc): 



F(z) = ^^oo(fc)^'. 



(12) 



k=0 



The de gree distribution for finite iV can be calculated as 
follows [241: 



p{k) dZ{N) 
NZ{N) dp{k) 



dz ^k- 

27Ti 



1-Nk T^N 



TToo 

* 

no^{k)w{N,k) 



(13) 



where the cutoff function w{N, k) is given by the above 
ratio of contour integrals. The integrals cannot be in 
general performed analytically for arbitrary A''. How- 
ever, the integration can be easily done numerically for 
different distributions iTooik) and sizes N as follows. We 



choose the contour of integration to be z 



re 



i4> 



with 



(j) £ (— TT, tt) and r smaller than the radius of convergence 
of F{z). Because the integrated function becomes con- 
centrated around zero for large iV, one does not need 
to integrate over the whole range of (j). We have writ- 
ten a procedure in Mathematica, which finds the value 
of r, for which the integrated function has the broadest 
maximum at (p = 0. This speeds up the convergence of 
numerical integration. Then, the procedure searches for 
0max for which the function falls to a sufficiently small 
part (typically 10~^°) of its maximal value. Then the 
function is integrated by means of the adaptive method 
over the range (— 0max, 0max)- Let us consider first the 
distribution ([5]) for 7 = 8, when it reduces to 



.(fc) 



A:(fc + l)(fc + 2)' 



for which 



Fiz) = 



z{3z - 2) - 2(z - l)2ln(l - z) 



(14) 



(15) 



so that the radius of convergence is one. Using the above 
numerical procedure we have calculated w{N, k) for iV = 
100, 200, 51200, shown in Fig. [iji. The theoretical 
value of a, predicted for this model, should be 1/2 [l7| . 
the same value comes from the correspondence to the 
zero-range process [2^. Therefore, for large TV, plots of 
w{N, k) for different sizes should collapse into a single 
curve m the rescaled variable x = k/N" = k/N^/"^. To 
check how fast the scaling is approached when N grows. 
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FIG. 1: (a) Plots of w{N, k) for multigraphs of sizes between 
100 (the narrowest peak) and 51200 (the widest peak), (b) 
The same data but plotted against the rescaled variable x = 
k/^NlnNf/^ . The highest peak shows data for iV = 100. 



instead of using the cutoff function directly it is better 
to consider its moments which can be calculated as 
follows: 



where the auxiliary function fm{z) is defined as 

d\™ 1 



frn{z) 



dz 



1- z 



(16) 



(17) 



We evaluated numerically the integrals in Eq. I|T6|) us- 
ing the same procedure as for the cutoff function, for 
m = 1,2,3 and sizes N = (1,2, . . . ,2^^) x 1000. The 
last value of TV is « 5 x 10^. In Fig. [2] we show plots 
of fj,rn{N)/N'-"^^^'>^^, where the sub-leading behavior has 
been exposed by dividing the data by the asymptotic 
form of iim 7V"("*+i). Since A'' spans several orders 
of magnitude, the horizontal axis has a logarithmic scale 
in order to make all points visible. The first moment 
^i{N)/N for the GNR model of growing tree network, 
calculated as in Ref. |[l^, is also plotted for compari- 
son. It approaches its maximal value as ^ N~^-^, which 
agrees with Eq. But the moments for the degen- 

erated graph show quite distinct behavior, namely they 
grow with IniV linearly or faster. This means that the 
subleading, power-like term from Eq. (flOl) might 

turn into a leading behavior which takes place for /3 — > 
and leads to a multiplicative, logarithmic correction. To 
check this, we fitted the formula 



fim{N)/N^'^+'^/^ = A,„{lnNf- + B.„ 



(18) 



to data points in Fig. [21 The fits are represented there as 
solid lines. The ratio of obtained values Ci : C2 ■ C3 — 
1.07 : 1.54 : 2.08 « 1 : 1.5 : 2 suggests that the moments 
behave as 



In other words, the cutoff scales as ~ {NlnNY^'^ . This 
is also confirmed in Fig. [TJd where the cutoff function is 
plotted in the rescaled variable x — k/{N\nNy^'^ . Ac- 
cording to our knowledge, this logarithmic correction has 
not been observed before. Its origin cannot lie in the dis- 
tribution 7roo(fc) only, because the growing networks with 
the same 7Too{k) do not have it, but it is the property of 
equilibrated graphs. Indeed, one can predict this scaling 
analytically, studying the cutoff function w{N, k): 



w{N,k) W{N,k)/W{N,Q), 



where 



W{N, k) 



dz 

2TTi' 



,k-l-Nk irN 



F-{z), 



(20) 



(21) 



with k = 2 and F{z) given by Eq. (fT5|) . Following the 
lines of Sec. 6 from Ref. [2^ one can argue that the func- 
tion under the integral is localized around z = 1 and 
can be expanded at this point. Choosing the contour of 
integration z — e*"^ we obtain 



W{N, k) 



rf0^W[iM+^2(3_„^1^^2)] ^^^^ 



with e small enough. Inserting now k = x{N hiNY/"^ 
and <j) = t{N\nNy^/'^ we have 

W{N,x{N\nNf'^)^{N\nN)-^/^ [ ^e't-t'+- 

J 27r 

(23) 

where higher terms are of order In In N/ In N and can be 
neglected for large TV. Then the above formula becomes 
a Gaussian integral and hence the cutoff function reads 



w{N, fc) « e 



* V (JV In JV)l/2 



(24) 



(19) 



with the predicted scaling. 

Let us now go to 7 > 3. In this region, the exponent 
a should be again 1/2 [2^. We calculated the first three 
moments ^m[N) for 7 — 3.5 from Eq. iflG]) for the same 
sizes iV as previously. Next, we fitted formula (flO)) to 
the data points. The best fit gives (3 ~ 0.24 for all three 
moments. This value is significantly greater than zero, 
thus the correction is now clearly of the power-like type. 
But the exponent /? « 1/4 is small enough to produce 
large corrections even for moderate-size networks. For 
example, if = IC*^, the correction is of order 0.16, with 
b being usually much larger than 1. But /3 grows with 
7, thus the convergence to the thermodynamic limit be- 
comes faster, e.g., for 7 = 4 we have estimated (3 > 0.4. 

We investigated also the range 2 < 7 < 3. Again, we 
calculated moments for different values of N for 7 = 2.5 
and fitted the formula (fTO| . The theory [23| predicts 
a — 1/(7—1) — 2/3 for this case. The agreement between 
the data and Eq. (fTO| is very good for ~ 0.33. 

At last, in order to check whether the formula Q for 
the large- a; behavior of the cutoff function holds for multi- 
graphs, we fitted a function A + ln(l + Bx) — (x/C)^ 
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with 



10 12 14 16 18 20 

IniV 
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FIG. 2: Plots of ^lrr,{N)/N'^'^+^^'^ calculated from Eq. ^ 
(squares), for degenerated graphs with 7 = 3 and for the 
GNR tree network with the same 7roo(fc) for comparison. Solid 
lines display functions: iflOll for lower-right panel, and ifTSjl for 
other panels, with parameters fitted to data points. 



to \n.w{x) obtained from numerical integration for N = 
100, 200, 51200, for 7 = 2.5, 3, 3.5. The assumed form 
of approximates the measured cutoff functions very 
well. The D's obtained for different sizes N tend to some 
limiting values which we found to be Dn-^oo = 2.0 for 
7 — 2.5, 1.9 for 7 = 3 and 3.1 for 7 = 3.5, with uncer- 
tainties of order 0.1. These values are in good agreement 
with the exponent 77 — 1/(1 — a) from Eq. ^ which gives 
2,2, and 3, respectively, and with theoretical results for 
the equivalent model [25l|. 

To summarize our findings for the degenerated graphs: 
for 7 = 3 we have observed logarithmic corrections to 
the scaling fccutoff ~ . The point 7 = 3 is the critical 
one; for 7 7^ 3, corrections are power-Hke, with the ex- 
ponent approaching zero for 7 ^ 3. The corrections are 
thus very strong for 7 ~ 3, even for very large networks. 
We have also examined the large-x behavior of the cut- 
off function w{x) and showed that it agrees with Eq. ^ 
derived for growing networks. This suggests some uni- 
versality, but we will see later that it does not hold for 
simple graphs. 



B. Equilibrated trees 

In Ref. the partition function Z{N) for equili- 

brated trees is found to be 



Z{N) = j> ^z-^-i^Gc(^). 



(25) 



Here Zqq{z) is a grand-canonical partition function 
obeying the equation 



F(z) = ^7r„o(fc + l)^ 



(27) 



k=0 



assuming that in the thermodynamic Hmit we want to 
obtain the degree distribution 7roo(fc). With help of 
Eq. (|26| . the partition function l(25|) can be rewritten 
as: 

Because we assume 7roo(0) — 0, the function F{z) equals 
\F[z), with F{z) as in Eq. ([12]) for multigraphs. This in 
turn leads to the following formula for Z{N): 

which is, up to a factor 1/iV, equivalent to the partition 
function (fTTj) for multigraphs with k = 2 — 1/N. This 
proves the equivalence between equilibrated multigraphs 
and trees, provided that the average degree is properly 
tuned, which has been already reported [22|. In our case, 
when L = N for multigraphs and therefore k = 2, a small 
difference in the values of moments fim arise. This does 
not, however, change the fact that in the thermodynamic 
limit both models become fully equivalent, and that for 
any finite TV the difference can be neglected. 

This means that the degree distribution for equili- 
brated trees behaves exactly as for multigraphs, that is 
we again have logarithmic corrections to the moments 
for 7 = 3, and power-law corrections for 7 7^ 3. 
This indicates also that the cutoff exponent a = 1/2 for 
7 > 3 and a = 1/(7 — 1) for 7 < 3. 



IV. MULTICANONICAL SIMULATIONS OF 
SIMPLE GRAPHS 

So far we have considered equilibrated multi- and 
tree graphs. But real-world networks are usually sim- 
ple graphs, that is they do not have multiple- and self 
connections, and have loops. Unfortunately, for these 
reasons simple graphs are not accessible with the tech- 
nique used before, because one does not know how to 
write the partition function as a single contour integral. 
Therefore, to obtain 7rAr(fc) one needs to turn to Monte 
Carlo (MC) techniques. Let us shortly describe here the 
general method which serves for this purpose, details can 
be found elsewhere [1, [2^ . We simulate graphs with fixed 
number of nodes N and links L. Each new graph gt+i is 
generated from the previous one gt by rewiring a single 
link. The move is accepted with the Metropolis proba- 
bility: 



Zgc{z) = zF{Zgc{z)), 



(26) 



Pigt 9t+i) = min <^ 1 



W{gt) 



(30) 
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which ensures that graphs are generated with correct 
weights W{g) from Eq. In order to restrict to simple 
graphs we reject moves introducing self- or multiple con- 
nections. This method allows for estimating 7r(fc) with 
good accuracy for N of order thousands, which is, how- 
ever, too small for our purpose. Larger networks are not 
accessible in a reasonable computer time. To show this, 
assume that 7r(fc) ~ k'"' and that we start from a network 
with all fc's much smaller than fccutoff ~ N". We assume 
that to obtain the experimental distribution TTcxpik) with 
acceptable level of noise, in the course of simulation each 
node has to change its degree between 1 and fccutoff many 
times. Let us denote by T the first-passage time from 
A; = 1 to fc = fccutoff. Our algorithm changes fc by ±1 
every step, so the process can be treated as a biased ran- 
dom walk in a potential V{k) = — ln7r(fc). Then T can 
be estimated as T '--^ ]\[<^i'y+'^) for 7 > 2, which is nothing 
more than the Arrhenius law T ^ exp [^(fccutoff )] with 
an additional correction factor. 

If, however, one could "flatten" 7r(fc) by appropriate 
reweighting of graphs, then instead of a biased random 
walk one would deal with a simple random walk where 
T ^ A^^". This is always better than the previous es- 
timate for the interesting range of 7. Moreover, by in- 
creasing the probability of graphs with fc > fccutoff, one 
could measure 7rcxp(fc) far above the cutoff. Here comes 
the idea of multicanonical simulations (MUCA) [27|, a 
similar idea has already been appHed to graphs [2a|. To 
apply MUCA in our case, we modify the weight function 
from Eq. ^ to: 

^(5)-^^np(fco, (31) 

where r(fci) is some function depending on the degree of 
one node, say the first one, and is chosen so that the dis- 
tribution 7r9at(fc), measured now only for the 1st node, 
is flat. The factor l/7roo(fci) gives additional advantage 
of flattening the distribution at no cost for k <^ fccutoff- 
To obtain the true distribution 7rcxp(fc) one multiplies 
7rffat(fc) by TToo{k)/r{k). Since all nodes are statistically 
equivalent in equihbrated graphs, this procedure must 
give the same result as standard MC simulations. The 
only question is to find the optimal r(fc). In our simu- 
lations, we applied the iterative method from Ref. [29| . 
which calculates r(fc) "on the fly" from sampled data. 
Because changing r(fc) during the simulation violates de- 
tailed balance and thus can change statistical weights, 
after obtaining sufficiently fiat histogram of 7rflat(fc) we 
fix r(fc) and perform a simple MC run with W{g) given 
by Eq. I|3ip . To speed up the simulation we divide the 
whole range of fc, for which we want to determine 7r(fc), 
into overlapping subranges of size « 50 and then glue re- 
sults. We vaHdated the method for graphs by comparing 
it with direct MC simulations for TV < 1000 (see Fig. [3]) 
and for multigraphs by comparing with exact results from 
Sec. III.A. 

To check how the cutoff scales with N , we simulated 




10 100 

k 



FIG. 3: Example of 7r(fe) for simple MC (squares) and MUCA 
(dotted line) for simple graphs with A'^ — 10^ and 7 = 3. Solid 
line shows iiaoik). Inset: the same for a broader range of k. 
MUCA simulation allows to estimate the tail of ttjv (fc) much 
better than simple MC. 



graphs with N = 625, 1250, 2500, . . . ,40000 nodes, and 
7 = 2.5, 3, 3.5, for the same noc{k) as before. In each 
case we performed 10^° to 10^^ MC steps (in units of 
single-link rewirings) for every subrange of fc, after fix- 
ing r{k). Each simulation was repeated 3-4 times in 
order to estimate statistical errors. We then calculated 
the moments fim of the cutoff function and checked if 
one can fit a simple power law fim — aN°'^™^-^^ to data 
points. Results are shown in Table Jl left, together with 
reduced values and significance levels ("Q" values) for 
the fits. We see that none of these cases can be accepted 
with the confidence level 0.05 which is commonly used 
in hypotheses testing. On the other hand, if we assume 
a more sophisticated form of fJ,m{N) with the sublead- 
ing correction from Eq. ifTOj) and, according to Ref. [lB |. 
with a(7 > 3) = 1/2, 0(7 = 2.5) = 1/(5 - 7) = 0.4, 
the agreement is very good, see Table HI right. We thus 
conclude that the subleading term plays an important 
role for moderate sizes. In Fig. [4] we show plots of ^ii{N) 
and both fitted functions. The subleading correction co- 
efficient b is very large. For instance, for 7 = 3.5 the 
correction is of order 0.1 for N ~ 10^, much larger than 
for multigraphs, perhaps due to much stronger structural 
constraints. Therefore, simple graphs approach the ther- 
modynamic limit even slower than equilibrated trees or 
degenerated graphs. 



Finally, we examined the large- a; behavior of w{x) us- 
ing the same method as for multigraphs. We found that 
the exponent 77 in Eq. |[9]) is not 1/(1 — a) as for multi- 
graphs but approximately 2.0 for all values of 7. This 
agrees with results of Ref. that the cutoff is always 
Gaussian in simple graphs, but it means also that t] is 
not universal. 
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FIG. 4: Plots of i^i{N) for simple graphs with 7 — 2.5, 3, 3.5. 
Solid and dotted lines are fits from Table U left and right, 
respectively. 



degenerated graphs simple graphs 




fcmax/^max ^max/^max 



FIG. 5: Left: rescaled plots of P(fcmax) for degenerated 
graphs, for 7 = 2.5 (top), 3 (middle) and 3.5 (bottom pic- 
ture), and for different sizes ranging from = 1000 up to 
32000 (1000: square, 2000: circle, 4000: triangle up, 8000: 

triangle down, 16000: diamond, 32000: triangle left). Solid 

—7+1 

lines are Prechet distributions P{x) oc x~'^e~^ fitted to 
data points for one chosen size A'' (usually A'' = 4000). Right: 
the same for simple graphs, lines are now: top and bottom 
- Gumbel, middle - Gumbel (thin line) and Frechet (thick 
line), fitted to black circles (iV = 2000). 



TABLE I: Results of fitting of the two different power laws 
(without and with a correction) to the measured moments Hm 
for simple graphs. For each 7, three consecutive rows present 
results for m — 1,2,3. "Parameters" indicate the set of free 
parameters during the fitting procedure, Q is the significance 
level of the fit. If a is marked as "fixed", its value is 0.4 for 
7 = 2.5 and 0.5 for the two other cases. 





parameters a, a 


=aiV"(™+i)(l + 6Ar-'') 
parameters a, b, /3, fixed a 


7 


a red. Q 


P h red. Q 


2.5 


0.366 4.5 <0.001 
0.373 3.5 0.004 
0.376 2.9 0.013 


0.39 5.4 2.40 0.05 
0.42 8 1.99 0.09 
0.44 11 1.72 0.14 


3 


0.392 2.64 0.021 
0.400 2.41 0.034 
0.405 2.34 0.039 


0.34 30 0.94 0.44 
0.40 80 0.87 0.48 
0.47 200 0.86 0.49 


3.5 


0.409 6.1 <0.001 
0.420 6.8 <0.001 
0.428 9.2 <0.001 


0.38 22 1.09 0.36 
0.48 56 0.73 0.57 
0.56 125 0.51 0.72 



V. CUTOFFS AND THE DISTRIBUTION OF 
MAXIMAL DEGREE 

So far we have studied the cutoffs in the degree distri- 
bution averaged over the ensemble of random graphs of 
a certain kind. The value A:cutoff ^ tells where the 
power law ends and TTN{k) starts to fall off rapidly. One 
can also pose the question: how does the maximal degree, 
^max, scale with N? Let a and a' be exponents in the 
power laws fccutoff ~ and fcmax ~ N" , respectively. 
The exponents a, a' may, but do not have to be equal, 
which has not been emphasized strongly enough in the 
literature on complex networks. To show this, let us con- 
sider the hypothetic network where all degrees ki,. . . ,kN 
are independent random variables taken from the dis- 
tribution p{k) ^ k^^ . The degree distribution 7r(A;) is 
then equal to p{k). In a real, finite network, degrees are 
(at least sHghtly) correlated due to structural constraints 
and 7r(fc) has a cutoff. In our simplified model, we can 
mimic these constraints by introducing the cutoff explic- 
itly: p{k) '-^ k^^w{^^) with some exponent a. Let us 
focus on a particular form of the cutoff: 

w{x) = exp(-ax''), (32) 

with parameters a, -q. This form of the cutoff function 
seems to be very general: in case of growing networks 
it is an asymptotic behavior for large fc, we also checked 
numerically that it holds for equilibrated graphs. For 
large N , the distribution P(fci„ax) of the maximal degree. 



^'(fcmax) « A^p(fcmax) / p{k)dk , (33) 
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takes its maximal value at k^^^ being a solution to 

p{k)dk^-Np\k*^^,), (34) 
which for the assumed cutoff gives 

Ar«(inAr)i/'?, a< ^^^^ 

This means that, modulo the logarithmic correction, a' = 
min(a, 1/(7 — 1)). One can also show that -P(fcmax) for 
a > 1/(7 — 1) is given by the Frechet extreme value dis- 

— '7 + 1 

tribution: P{x) oc x~'^e~'^ , while it approaches the 
Gumbel distribution: P{x) oc e~^~'^ for a < 1/(7—1), 
in the properly rescaled variable x = A + Bfcmax- Be- 
cause these two distributions are narrow, fc^ax approx- 
imates also the mean value (/cmax)- We can learn two 
things from Eq. l(35|) . First, fcmax scales as a pure power 
of TV with the same exponent as for fccutoff, a' — a, only 
if a = 1/(7 — 1). The distribution of /cmax is then neither 
Frechet nor Gumbel, but has a more complicated form, 
as follows from Eq. ((33| : 



P{x) = Mx ''w{x) exp 



M y ^w{y)dy 



(36) 



where A/" is a normalization coefficient and x = 
kniax/N^/'-''~^\ and depends on the exact form of w{x). 
Second, a' may be smaller than a, that is the maximal 
degree grows slower than the cutoff. At first sight this 
may appear counter-intuitive: one could have expected 
that 7r(fc) for large k is dominated by the maximal de- 
gree distribution and thus kmax must not grow slower 
than fccutoff- 

The formula for fc^ax presented above works surpris- 
ingly well for real graphs where the cutoff stems from 
correlations between nodes degrees. In order to check 
Eq. l(35|) . we performed multicanonical simulations and 
obtained distributions /'(fcmax) for multigraphs and sim- 
ple graphs with 7 = 2.5, 3 and 3.5. In light of what has 
been said in Sec. III.B, the scaling of fcmax for trees has to 
be identical to that for multigraphs, so it is not necessary 
to consider trees separately. The procedure was similar 
to that described in the previous section, with the only 
difference that we fiatten P(fcmax) and not 7r(fc). In each 
simulation we performed around 10* Monte Carlo steps 
after fixing the weight r(fci„ax)- 

Let us discuss first multigraphs. In Table [ill left, 
we present values of a, ry for various 7, which we con- 
firmed numerically in previous sections. The fourth col- 
umn shows scaling laws for fej^^x predicted by means of 
Eq. l(35|l . In Fig. [5] we plotted the experimental distribu- 
tions P(A;i„ax), in the rescaled variable x = fcrnax/^max- 
Plots for different N show good agreement of positions 
of the maximum. The distribution for 7 = 3.5 is also 
approximately Frechet, which agrees with recent find- 
ings |33|, the two for 7 = 2.5, 3 deviate slightly from 



TABLE II: Values of exponents a, rj, and scaling of fc^ax 
for degenerated and simple graphs, for different power laws 
■K{k) ~ k Assumed scalings are: fccutoff ~ N" , \nw{x) ~ 
—x^. The exponents a,rj have been confirmed numerically 
in the previous section. The formulas for fc^ax are obtained 
from Eq. 
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multigraphs 


simple graphs 




^ V fcmax 


a 


V ^max 


(2,3) 


7-1 7-2 


1 

5-7 


2 iVi/(5-^'(lniV)i/2 


3 


1/2+log 2 


1/2 


2 iV^/^ 


(3,00) 


1/2 2 iVi/(^-i) 


1/2 


2 iVi/{7-i) 



Frechet in the tail. The perfect scaling of fcmax means 
that the thermodynamic limit for the distribution of max- 
imal degree is approached much faster than for 7r(fc), and 
is reached already for N ^ 10^ nodes; no subleading cor- 
rections are necessary. 

We repeated the same procedure for graphs, see Table 
HU right, and Fig. \E[ Again, one sees very good agree- 
ment, except for 7 = 3, where one obtains better collapse 
for a sHghtly different value fcj^ax N'^''^^- The shape of 
-P(fcinax) is better approximated by Gumbel than Frechet 
distribution. Also for 7 = 3.5 and N < 32000, the dis- 
tribution -P(fcmax) deviates from Frechet and is closer to 
Gumbel. 

To conclude, we have shown that the distribution of 
the maximal degree exhibits the scaling fcmax ~ N°' but 
in general with a different exponent than that in the cut- 
off. The large-iV scaHng is, however, approached much 
faster than in case of the degree distribution and sublead- 
ing corrections can be safely neglected in most cases for 
networks of order IC* nodes. It is also worth mentioning 
that a' depends on both the exponents 7 and a. 



VI. SUMMARY 

In this paper we considered finite-size effects in the de- 
gree distribution of equilibrated networks. We showed 
that the convergence towards the thermodynamic limit 
is very slow, thus in order to get reasonable estimation of 
the cutoff one has to consider subleading corrections to 
the scaling fccutoff ~ N"". For multigraphs and 7 = 3, the 
correction turns into a leading behavior ~ (TV In TV) 
and therefore does not vanish for ^ 00. This has 
to be taken into account when comparing any numeri- 
cal results for finite networks to that derived analytically 
for the leading behavior only. We checked also that the 
asymptotic behavior of the cutoff function is very sim- 
ple: hiw{x) ^ —x^, but r] is not universal among differ- 
ent classes of graphs. We argued that the maximal de- 
gree fcmax reaches the asymptotic power-law scaling faster 
than the cutoff but scales in many cases differently than 

fccutoff • 

The networks presented here do not have explicit 
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degree-degree correlations. It would be interesting to 
check which of the results will survive the introduction 
of correlations, i.e., are universal to some extent. 
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